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ABSTRACT 

We launch a fully relativistic study of the formation of supermassive black holes via the collapse of 
supermassive stars. Here we initiate our investigation by analyzing the secular evolution of supermassive 
stars up to the onset of dynamical instability and collapse. We focus on the effects of rotation, assumed 
uniform, and general relativity. We identify the critical configuration at which radial instability sets in 
and determine its structure in detail. We show that the key nondimensional ratios R/M, T/\W\ and 
J/M 2 (T is the rotational kinetic energy and W is the gravitational potential energy) for this critical 
configuration are universal numbers, independent of the mass, spin, radius or history of the star. We 
compare results from an approximate, analytic treatment with a fully relativistic, numerical calculation 
and find good agreement. We solve analytically for the time evolution of these parameters up to the onset 
of instability. Cooling by photon radiation drives the evolution, which is accompanied by mass, angular 
momentum and entropy loss. The critical configuration serves as initial data for a future relativistic, 
hydrodynamical, 3D simulation of the collapse of an unstable supermassive star. Since this implosion 
starts from a universal critical configuration, the collapse is also uniquely determined and should produce 
a universal gravitational waveform. In this paper we briefly speculate on the possible outcome of this 
collapse and asses to what extent it offers a promising route to forming a supermassive black hole. 



1. INTRODUCTION 

Recent observations provide increasingly strong evi- 
dence that supermassive black holes (SMBHs) exist and 
that they are the sources that power active galactic nuclei 
(AGNs) and quasars (see, e.g., Rees 1998, for a summary 
and references). For example, the 1.3 cm maser-emission 
line of H2O allows for an extremely accurate mapping of 
gas motions in a disk at the core of the galaxy NGC 4258. 
The inner edge of the disk has a radius of ~ 0.1 pc, and 
the inferred velocity profile is consistent with Keplerian 
motion around a central object of 3.6 x 10 7 M Q (Claussen 
& Lo 1986; Watson & Wallin 1994; Miyoshi et. al. 1995). 
The most conservative candidate for so massive an object 
in so small a volume is a black hole. 

However, the scenario by which SMBHs form with 
masses in the range of 1O 6 -1O 1O M0 is still uncertain (see 
Rees 1984 for an overview). Viable stellar dynamical and a 
hydrodynamical routes leading to the formation of SMBHs 
have both been proposed. In one stellar dynamical sce- 
nario, the primordial gas first fragments into stars, which 
form a dense stellar cluster or galaxy core. Repeated col- 
lisions and mergers of the stars in the cluster lead to the 
build-up of massive stars, and the subsequent collapse of 
the massive stars result in the formation of stellar-mass 
black holes. As these holes continue to swallow stars and 
grow, they settle toward the galaxy center, merge and trig- 
ger the build-up of one or more SMBHs (see, e.g., Quinlan 
& Shapiro 1990 and references therein). In another dy- 
namical scenario, the gravothermal catastrophe (secular 
core collapse) in a dense cluster composed of neutron stars 
or stellar-mass black holes may drive the cluster core to a 
relativistic state. Sufficiently relativistic clusters are dy- 
namically unstable to catastrophic collapse and may give 



birth to SMBHs (Zel'dovich & Podurets 1965; Shapiro & 
Teukolsky 1985a, 1985b; Quinlan & Shapiro 1987). It may 
be more likely that the contracting primordial gas builds 
up sufficient radiation pressure to inhibit further fragmen- 
tation and hence prevent the formation of a stellar clus- 
ter. Alternatively, a supermassive gas cloud may build 
up from the fragmentation of stars in collisions in stel- 
lar clusters (Sanders 1970; Begelman & Rees 1978). In 
such hydrodynamical scenarios, a massive black hole may 
form directly from gas, perhaps after a phase in which the 
gas evolves in a quasistationary manner as a supermassive 
star (SMS). Such a phase can only be transient, since ulti- 
mately general relativity will drive such stars dynamically 
unstable (Iben 1963; Chandrasekhar 1964a, b; Feynman, 
unpublished, as quoted in Fowler 1964). The structure 
and evolution of supermassive stars up to the onset of in- 
stability has been explored previously by many authors, 
including Hoyle & Fowler (1963a, b), Fowler (1964, 1966), 
Bisnovatyi-Kogan, Zel'dovich & Novikov (1967), Appen- 
zeller & Fricke (1972) and Fuller, Woosley and Weaver 
(1986) (see Zel'dovich & Novikov, 1971, and Shapiro and 
Teukolsky, 1983, for reviews and further references). 

In addition to their importance for a fundamental un- 
derstanding of AGNs and quasars, supermassive objects 
and the formation of SMBHs have recently attracted new 
interest because they are likely candidates for detection by 
currently proposed space-based gravitational wave detec- 
tors like the Laser Interferometer Space Antenna (LISA). 
Because of its long arm length, this detector would be very 
sensitive to long wavelength and low frequency radiation, 
and therefore supermassive objects are among the most 
promising sources (see, e.g., Thorne & Braginsky 1976; 
Thorne 1995). In particular, LISA may be able to detect 
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the collapse of a SMS to a SMBH. Even more promising 
is the possible detection of the coalescence of two SMBHs 
(LISA Pre-Phase A report 1995). The likelihood of such 
an event, however, largely depends on how SMBHs form 
and is therefore still uncertain. 

In a series of papers we revisit the formation of SMBHs 
via the collapse of SMSs, focussing on the influence of ro- 
tation and general relativity. We analyze the secular con- 
traction of a uniformly rotating equilibrium configuration 
via thermal emission and mass loss. We concentrate on 
a configuration rotating at the mass-shedding limit. In 
Baumgarte & Shapiro (1999, hereafter Paper I), we have 
shown that the luminosity from such a star is considerably 
reduced below the value of a nonrotating spherical star of 
the same mass. 

In this paper, we analyze the structure and stability 
against collapse of fully relativistic, rotating n = 3 poly- 
tropes in stationary equilibrium. SMSs to which these cal- 
culations apply are radiation-dominated, isentropic con- 
figurations of sufficient mass that neither nuclear burning 
nor electron-positron pairs are important before the stars 
reaches the onset of relativistic gravitational instability. 
Stars with M ;> 10 6 M Q fall in this category (Zel'dovich & 
Novikov 1971; Fuller, Woosley & Weaver 1986). Moreover, 
the evolutionary timescale due to cooling has to be longer 
than the hydrodynamic timescale for the star to evolve in 
a quasistationary fashion. According to equations (9) and 
(10) below, we find that this constraint is satisfied for all 
masses M <; 1O 13 M . 

We track the quasistationary evolution of such stars to 
the point of onset of radial instability. We first identify 
this critical configuration analytically by means of a post- 
Newtonian energy variational method. As pointed out 
by Bisnovatyi-Kogan, Zel'dovich & Novikov (1967), de- 
termining the onset of relativistic collapse in a rapidly ro- 
tating SMS requires a second-order post-Newtonian anal- 
ysis, not a first-order, as in the case of a nonrotating 
star. We provide an approximate, post-Newtonian ana- 
lytic analysis (rigorizing and completing the qualitative ar- 
gument of Bisnovatyi-Kogan, Zel'dovich & Novikov (1967) 
and Zel'dovich & Novikov (1971), who only provide a di- 
mensional estimate of the second post- Newtonian term). 
Then, to establish the result rigorously, we use a recent 
numerical code (Cook, Shapiro & Teukolsky 1992, 1994) 
to construct rotating equilibrium models in full general 
relativity. We show that the mass of the critical config- 
uration is a unique function of the specific entropy. The 
values of R/M, T/\W\ and J/M 2 at the onset of collapse 
are universal numbers, independent of the mass or prior 
evolution. Here M is the mass, R the (polar) radius, J the 
angular momentum, T the rotational and W the gravita- 
tional potential energy of the star. We also speculate on 
the likely outcome of collapse for stars which do not dis- 
rupt due to thermonuclear explosions during collapse; it is 
found that stars with M > 10 5 M e and initial metallici- 
ties Z < 0.005 do not explode (Fuller, Woosley & Weaver 
1986). Since these stars start collapsing from a univer- 
sal critical configuration, the subsequent collapse is also 
uniquely determined and should produce a unique gravi- 
tational waveform. We postpone a detailed discussion of 
this phase for a future paper in which we follow the dynam- 
ical collapse numerically in general relativity (Baumgarte, 



Shapiro & Shibata 1999). 

The key goals of our study are to decide whether a 
SMBH can really emerge from the collapse of a SMS and to 
determine the hole parameters if indeed it can be formed 
this way. Alternatively, a rotating supermassive cloud or 
star could collapse to a weakly relativistic disk (e.g. Wag- 
oner 1969; Loeb & Rasio 1994). If the collapsing innermost 
region enters the strong-field domain, the angular momen- 
tum of this matter must be below the maximum value of 
a Kerr hole (J/M 2 = 1) for black hole formation to occur 
eventually. What happens if the angular momentum ex- 
ceeds this limit? Does angular momentum dissipation by 
outflowing gas allow for black hole formation of the core? 
Or, does gravitational radiation, following the formation 
of bars or axial currents carry away enough angular mo- 
mentum to permit collapse? Finally, if a SMBH can form 
by the collapse of a SMS, what is its mass and spin given 
the mass and spin of the SMS at the onset of collapse? 

In this paper we deal primarily with the structure, sta- 
bility and early secular evolution phases of the SMS sce- 
nario. Our calculation, in effect, sets up the initial data at 
the onset of collapse. Tracking the subsequent dynamical 
evolution of these initial data will resolve the key ques- 
tions posed above, and in this paper we will only speculate 
briefly on the outcome of the dynamical collapse. 

This paper is organized as follows: In Section ^ we pro- 
vide a qualitative overview of the problem and present our 
basic assumptions. In Section || we discuss the equilibrium 
and stability of rotating, relativistic SMSs. In particular, 
we determine the critical configuration at which an evolv- 
ing SMS becomes dynamically unstable to radial pertur- 
bations. We compare re sult s from an approximate analyt- 
ical treatment (Section |3.2| ) with tho se f rom a numerical, 
fully relativistic calculation (Section |3.3| ) . Having identi- 
fied the onset of instability, we then solve analytically for 
the evolution of SMSs during the secular contraction phase 
up to this critical configuration in Section || In Section || 
we provide some qualitative arguments which suggest that 
the direct formation of SMBHs from the collapse of SMSs 
indeed may be possible. We summarize and discuss our 
results in Section ^. Except where noted otherwise, we 
adopt geometrized units with c = 1 = G throughout this 
paper. 

2. QUALITATIVE OVERVIEW AND BASIC ASSUMPTIONS 

SMSs may form if collapsing primordial gas builds up 
enough entropy so that the radiation pressure can slow 
down the collapse (see Begelman & Rees, 1978, for an al- 
ternative scenario). Further contraction will then spin up 
the newly formed SMS to the mass-shedding limit, pro- 
vided that the gas had some initial angular momentum and 
that viscosity maintains uniform rotation. The SMS will 
then evolve secularly along the mass-shedding limit, simul- 
taneously emitting radiation, matter and angular momen- 
tum (see, e.g., Bisnovatyi-Kogan, Zel'dovich k. Novikov 
1967; Zel'dovich & Novikov 1971). Once it reaches the on- 
set of radial instability, the star collapses on a dynamical 
timescale, and may ultimately form a SMBH. 

For sufficiently massive objects (M ^ 10 6 Af Q ), the 
equation of state is dominated by thermal radiation pres- 
sure. It can also be shown that SMSs are convective (see 
Loeb & Rasio 1994, for a simple proof) with constant en- 
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tropy per baryon, 
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3 n b 
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where nb is the baryon density and a is the radiation den- 
sity constant. These conditions imply that the structure 
of a SMS is that of an n = 3 polytrope 



P = Kp 4 ' 3 , 



where 



K = K{s) = 



a ( 3s 



3 \ 4m#a 



4/3 



= const. 



(2) 



(3) 



(see eq. 17.2.6 in Shapiro & Teukolsky 1983). Here, mn 
is the mass of a hydrogen atom and K has been evaluated 
for a composition of pure ionized hydrogen. In Newtonian 
theory, the mass of a static, equilibrium n = 3 polytrope 
is uniquely determined by the polytropic constant K alone 



M = K 



3/2 



(4) 



(the numerical coefficients k\ and ki, derived from Lane- 
Emden functions, are given in Table 1 below). Inverting 
this relation, we can write K in terms of M as 



# = 1.01 x 10 3 cm 2 / 3 f — 



2/3 



(5) 



Since n — 3 polytropes are extremely centrally condensed, 
the self-gravity of the outermost envelope of a uniformly 
rotating configuration ca n be neglected and the Roche ap- 
proximation (see Section 3.1) can be applied. 

In Newtonian gravitation, n = 3 polytropes are margin- 
ally stable to radial collapse. However, even very small 
general relativistic corrections make these polytropes un- 
stable, and some mechanism has to be invoked to prevent 
collapse. In this paper, we focus on rotation, which up to 
a critical configuration along the mass-shedding sequence 
can stabilize SMSs. In Section || we will identify this crit- 
ical configuration (which we will call "configuration A"). 
Alternatively, gas pressure may stabilize SMSs in the ab- 
sence of angular momentum (see Zel'dovich & Novikov 
1971; Shapiro & Teukolsky 1983), but even a small degree 
of rotation will dominate gas pressure (cf. Section 3.2). 
A dark matter background also tends to have a stabi- 
lizing effect (Fuller, Woosley & Weaver, 1986), but it is 
also much less important than rotation (Bisnovatyi-Kogan 
1998). Since the ratio between the rotational kinetic and 
the potential energy, T/|W|, is always very small for a 
uniformly rotating star during secular evolution along the 
mass-shedding sequence (see eq. JjiJ and Figure ||), it is 
unlikely that nonradial modes of instability are important 
during this phase. Once the star has become unstable and 
collapses, however, T/\W | can become very large, and it is 
likely that such nonradial modes develop (see Section ||) . 

Massive enough stars [M J> 10 6 M Q ) do not reach suf- 
ficiently high temperatures for nuclear burning to become 
imp ortan t before reaching the onset of instability (see Sec- 
tion 3.4). Also, electron-positron pairs play a negligible 
role in this regime (Zel'dovich & Novikov 1971). The 



evolution of the SMS along the mass-shedding sequence 
is then determined solely by cooling via thermal photon 
emission and loss of mass and angular momentum. As 
we have shown in Paper I, the luminosity of a rotating 
supermassive star at mass-shedding is 



L = 0.639 L Edd , 
where the Eddington luminosity is 

4ttM 



-^Edd = 



(6) 



(7) 



The opacity k is dominated by Thomson scattering off free 
electrons 



k t = 0.2(1 + A 



H 



cm 2 g 



(8) 



where Xh is the hydrogen mass fraction. For M J> 
10 5 Af Q , temperatures are low enough that Klein-Nishina 
corrections can be neglected (see, e.g., Fuller, Woosley & 
Weaver 1986). 

As we will find in Section |], the evolutionary timescale 
for secular evolution along the mass-shedding sequence is 
given by 



W = 8.8x10"b(^)' 



(9) 



where the coefficient has been evaluated near the critical 



configuration (see eq. (67)). The dynamical timescale is 



*dyn ~ (G Pc )- 1/2 



(10) 



1.3 x 10" 2 s 



M 



7 x 10- 9 



-1/2 



where we have used eq. ( p9| ) below to evaluate p c in terms 
of the non-dimensional density p c . Given the critical val- 
ues for p c and R/M that we find in Section |^, the evolu- 
tionary timescale is longer than the dynamical timescale 
for M < 6.7 x 10 13 M Q . This mass is larger than the value 
derived in Shapiro & Teukolsky (1983), where a similar 
estimate has been made for nonrotating SMSs stabilized 
by gas pressure. 

Also important is the viscous timescale, which is very 
uncertain. Not surprisingly, the microscopic viscosity due 
either to collisions between ions or to radiation yields 
timcscales which are larger than the evolutionary timescale 
by many orders of magnitude (see Kippenhahn & Weigert 
1990). The effect of turbulent viscosity can be estimated 
by assuming that the velocity of the turbulent motion Vt 
is an appreciable fraction of the velocity of sound, 



V t = abound, 



(11) 



where we take the dimensionless viscosity parameter a to 
lie in the range 

0.01 <a<l (12) 

(see Shakura & Sunyaev 1973; also Zel'dovich & Novikov 
1971; Balbus & Hawley 1991). Assuming the characteris- 
tic scale of nonuniformity to be some fraction of the stellar 
radius i?, we can then estimate the turbulent viscosity to 
be 

n t ~ pRv t = pRav sound . (13) 
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The corresponding timescale is then 
R 2 p R 

°f]t abound 



a 1 td y 



(14) 



It has been argued that both convection and magnetic 
fields serve to generate such a turbulent viscosity in SMSs 
(Bisnovatyi-Kogan, Zel'dovich & Novikov 1967; Wagoner 
1969). 

For SMSs with masses M < 6.7 x 1Q U M , we therefore 
find the hierarchy 



^dyn ^vis < ^ ^cvol • 



(15) 



Eq. (15) justifies our assumption that the star evolves 
in a quasistationary manner along the mass-shedding se- 
quence. The last inequality implies that turbulent viscos- 
ity, if present, will keep the SMS in uniform rotation. We 
will assume this result in treating the secular evolution up 
to the onset of instability. Once the star becomes unstable 
and starts collapsing, however, the evolution will proceed 
on the dynamical timescale, and it is unlikely that the 
viscosity will be efficient enough to maintain uniform ro- 
tation. It is more likely that cylindrical mass shells will 
conserve their angular momentum, which we assume for 
our qualitative arguments in Section [s]. 

3. EQUILIBRIUM AND STABILITY 

In this Section we analyze the equilibrium and stability 
of rotating SMSs and identify the critical configuration at 
which radial instability sets in. We briefly revi ew the pre- 
dictions of the Roche approximation in Section ft.l| . In Sec- 
tion 3.2 we present an analytic, post-Newtonian model cal- 



culation, which provides qualitative insight into the scaling 
behavior and approximate num eric al parameters for the 
critical configuration. In Section |3.3| , we present numerical 
results from a fully relativistic calculation, which allows us 
to identify the critical configuration and its characteristic 
parameters (including J/M 2 and R/M) independently of 
the approximations used in the analytical treatment (in- 
cluding the post-Newtonian and Roche approximation). 
We find remarkable agr eem ent between the two calcula- 
tions. Lastly, in Section 3.4, we restore cgs units and pro- 
vide expressions for the physical parameters in SMSs. 

3.1. Predictions of the Roche Model 

The Roche approximation for rotating n = 3 polytropes 
has been derived and discussed in the literature (see, e.g., 
Zel'dovich & Novikov 1971; Shapiro & Teukolsky 1983). 
Here we briefly summarize the key elements for applica- 
tion below (see also Paper I). 

Stars with soft equations of state are extremely centrally 
condensed: they have an extended, low density envelope, 
with the bulk of the mass concentrated in the core. For 
ann = 3 polytrope, for example, the ratio between central 
density to average density is p c /(p) — 54.2. The gravita- 
tional forces in the envelope are therefore dominated by 
the massive core, and it is appropriate to neglect the self- 
gravity of the envelope. In this case, it is easy to show 
that matter at the equator of a rotating star orbits with 
the Kepler frequency when the ratio between the equato- 
rial radius R e and the polar radius R p reaches 



P / shcdd 



(16) 



At this point, the angular velocity of the star assumes its 
maximum value before mass-shedding from the equator 
sets in: 

/2\ 3/2 / M \ 1/2 



The polar radius hardly changes from the value of the ra- 
dius of the nonrotating, spherical star, a result which has 
been verified numerically (Papaloizou & Whelan 1973). 
Since the bulk of the matter is concentrated at the core 
and hardly affected by the rotation, the moment of inertia 
of inertia of the star barely changes with rotation and is 
well approximated by the nonrotating value 



1= ~(r 2 )M = ^0.1130 R 2 p M, 



(18) 



where the brackets denote a mass weighted average over 
the star, and where the last equality holds for n — 3 
(this app roximation turns out to be less accurate; cf. Sec- 
tion [T^). The ratio between the kinetic and potential 
energy at mass-shedding then becomes 



T 

W\ 



shcdd 



(V2) ^ "shcdd 

(3/2)M 2 /i? p 

2_Y 0.1130 

3 



0.00744. 



(19) 



This result predicts that T/|W| of a maximally rotating, 
n = 3 polytrope (or any star for which the Roche approx- 
imation applies) is a universal constant, i nde pendent of 
mass, radius, or orbital velocity. In Section 3.3 we will see 



that precise numerical models predict a somewhat larger 
value of T/|W|, which nevertheless also remains approxi- 
mately constant along a mass-shedding trajectory (see Fig- 
ure [I| below). 

3.2. Analytical Model 

To determine the equilibrium and stability of a rotat- 
ing SMS, we write its total energy as the sum of the in- 
ternal energy U, the potential energy W, the rotational 
energy T, a post-Newtonian correction Epn and a post- 
post- Newtonian correction Epp^- For a, n — 3 polytrope, 
these terms can be written 



E = kiKMp x J 3 - k 2 M 5/3 p 1 J 3 + 
k 3 j 2 M 7 / 3 p 2 J 3 - k 4 M 7 / 3 p 2 J 3 



(20) 



k 5 M 3 p c 



where we have defined j = J/M 2 and have neglected 
corrections due to deviations from sphericity. This ne- 
glect is justified, since these corrections scale with T/|W|, 
which according to ( |l9| ) is always very small. Even though 
the value of the post-post-Newtonian correction Eppn is 
very small, this term is crucial for determining the crit- 
ical, marginally stable configuration, as emphasized by 
Zel'dovich and Novikov (1971). The values of the nondi- 
mensional coefficients fcj are listed in Table 1; they are 
constructed from Lane-Emden functions. 

Note that for any polytrope K n l 2 has units of length. 
We can therefore introduce nondimensional coordinates by 



BAUMGARTE & SHAPIRO 



5 



Table 1 

Values of the structure coefficients for n = 3. 



fcl 


fc 2 


^3 


fc 4 


fc 5 


1.7558 1 


0.63899 1 


1.2041 1 


0.918294 2 


0.331211 3 



References. — (1) Lai, Rasio & Shapiro, 1993; (2) 
Shapiro & Teukolsky, 1983; (3) Lombardi, 1997 



setting K — 1 (see also Cook, Shapiro & Teukolsky 1992). 
We will denote values of nondimensional variables in these 
coordinates with a bar (for example M). Values of these 
quantities for any other value of K can be recovered easily 
by rescaling with an appropriate power of K n / 2 = K 3 / 2 ; 
for example M — K 3 I 2 M and p = K~ 3 p. 

Taking the first derivative of (12^) with respect to the 
central density yields a condition for equilibrium 



= — = fciM - k 2 M 5/3 + 2k 3 j 2 M 7/3 x 
ox 



2k A M 7/3 x - 3k 5 M 3 x 2 , 



(21) 



_l/3 

where x = p c ■ For stable equilibrium, the second deriva- 
tive of eq. ( pp|) has to be positive. A root of the second 
derivative therefore marks the onset of radial instability 







d 2 E 
dx 2 



2k 3 j 2 M 7/3 - 2k 4 M 7/3 - 6k 5 M 3 x. (22) 



Solving eq. (|2^) for x immediately yields 

k 3 j 2 - fe 4 



pT = - = 



3fc 5 M 2 / 3 



For stability, we obviously require p c > or 

1/2 



j > Jmin = ( ^ ) = °' 8733 - 



(23) 



(24) 



Note that stable configurations with j — j m j n have p c = 0. 

SMSs with angular momenta less than j m ; n can also be 
stabilized by the nonvanishing contribution of the plasma, 
which effectively decreases the polytropic index n to a 
value slightly smaller than 3 (see Shapiro & Teukolsky 
1983, Chapt. 17). Here, because we are interested in the 
effects of rotation, we neglect the plasma and approxi- 
mate supermassive stars as strict n = 3 radiation domi- 
nated polytropes. Our neglect of gas pressure breaks down 
for angular momenta less than j m i n , since it would make 
wrong predictions for the stability of supermassive stars. 
We expect that most equilibria which form in nature have 

j > jmin. 

Inserting (^3|) into (|2l]) , we find the mass of equilibrium 
stars at the onset of instability is 



M 2/3 = fci 



1 

3fc 5 



for 



(25) 



For j = jmin, and hence p c = 0, the configuration has infi- 
nite radius and we obviously recover the (unique) Newto- 
nian mass of a n = 3 polytrope 



Mmin = I 4 



3/2 



= 4.5548 



(26) 



(where we keep enough digits for purposes of comparison 
later on). This configuration corresponds to point B in 
Figs. U and 0. Note that from equation (^), all criti- 
cal configurations at the onset of instability have a mass 
greater or equal to M m in- 

The maximum angular velocity of a rotating star oc- 
curs at mass-shedding, where we know the value of T/|W| 
(eq. [l9] in the Roche approximation) . It is convenient to 
express T/|W| in terms of x, M and j 



T 

W\ 



k 3 j 



'■M 7 / 3 x 2 



k 3 j 2 M 2 / 3 i 



k 2 M 5 / 3 x 



(27) 



The last stable configuration for a given T/|1V| can now 
be found by substituting ( p3| ) into (|27|), which yields a 
quadratic equation for j 2 : 



k/\ o fc^fcs T 

Jin 2 _ o. a — n 

fc 3 J fc 2 \w\- u ' 




fc 2 \wy 



(28) 



(29) 



(The second solution can be disregarded, since it gives a 
negative value for j 2 ). To first order in T/|1T|, this result 
can be approximated by 



. _ . h 3 fc 2 fc 5 T 



j min (l + 0.377^). (30) 



The mass of this critical configuration now can be found 
by inserting (pp|) into (25), which yields 



M 2 ' 3 



fci 
k 2 



fci w 



(31) 
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or, to first order in (T/|VT 



M — M [l + 9 -^l(JL^ 



Mr, 



t T \ 2 



(32) 



Adopting the Roche prediction (|l9|) for T/|W|, we find 



j crit = 0.8757 
M cri t =4.5551. 

Inserting these values into (^3|) yields 
Oc)crit = 6.6 x 10-°. 



(33) 
(34) 



(35) 



These values correspond to the critical configuration A in 
Figs. | and | 

It is obvious from expressions (|3Cl| ) and (|3^) that the 
critical values j cr it and M C rit depend on several approxi- 
mations. In this section, we have adopted the Roche pre- 
diction © for T/\W\ at mass-shedding, even though self- 
consistent numerical calc ulati ons yield a value that differs 
by about 20% (cf. Section 3^3 and Table 2). Moreover, the 
critical values depend on the second post-Newtonian coef- 
ficient fcs. We have adopted a value that only takes into 
account second post-Newtonian corrections to the energy 
of a spherical, nonrotating configuration (Lombardi 1997), 
and have neglected first post-Newtonian corrections to the 
Newtonian rotational kinetic energy. The later would have 
the same form as the second post-Newtonian term Eppn, 
and would therefore introduce an unknown correction to 
the coefficient k$ (which we expect to be small for small 
values of T/|IF|). Together, this means that our analytic 
calculation can only yield approximate values for the crit- 
ical configuration. 

It is interesting to note, however, that our analytic 
model predicts j cr it to be less than unity, which is the max- 
imum value allowed for a Kerr black hole. A supermassive 
star that evolves along the mass-shedding sequence and 
ultimately becomes unstable could therefore collapse to a 
black hole without having to loose additional angular mo- 
mentum (see also Section ||). 

Lastly, we estimate the compaction of the critical con- 
figuration by rewriting equation (^7|) in terms of the polar 
radius R p . Neglecting deviations from sphericity, we find 



T 

W\ 



2J Rp 
6IM 2 



1 



M 



2 x 0.1130 i? K 



or 



Rp 
Jl 



= 4.425 



T/\W\ 



(36) 



(37) 



where we have dropped the bars on the left hand side, 
since R p /M is a dimensionless combination, and have used 
eq. (|l8|) to evaluate /. Inserting (K^j and (p3|), we find 



Rp 
~M 



456. 



(38) 



This shows that the critical configuration is only very 
mildly relativistic, and that a post-Newtonian approxima- 
tion is adequate for determining its equilibrium structure. 



The small admixture of thermal gas pressure can sta- 
bilize a nonrotating SMS for values of R/M exceeding 
(i?/M)s as = 1.59 (M/Mq) 1 / 2 (eq. (17.4.11) in Shapiro & 
Teukolsky, 1983). For masses M £ 1O 5 M , this ratio is 
larger than that given eq. (|38|). Consequently, rotation 
serves to stabilize these masses beyond the point up to 
which thermal pressure alone can stabilize them. Hence 
thermal gas pressure can be ignored, as we have done, 
in determining the point of onset of instability in mass- 
shedding configurations with M ;> 1O 5 M . 

Note that all the nondimensional parameters of the crit- 
ical configuration, including T/|VT|, R/M and J/M 2 , are 
universal constants, and do not depend on the mass of the 
SMS or its history. 

3.3. Numerical Results 

In this Section we present numerical models of fully rel- 
ativistic, uniformly rotating n = 3 polytropes. We have 
constructed these models with the numerical code of Cook, 
Shapiro & Teukolsky (1992, 1994). As a typical result, we 
show density contours of the last stable (critical) configu- 
ration (configuration A) in a plane containing the axis of 
rotation in Figure [l]. 

In the static limit, we recover the Tolman-Oppenhei- 
mer-Volkoff result for relativistic hydrostatic equilibrium 
in spherical symmetry. We mark this curve by the long 
dashed line in Figures || and ||. Note that for a polytropic 
index n — 3, Newtonian gravity predicts a unique value 
for the mass M, which is independent of the central den- 
sity (equation (J2q)). In general relativity, M is no longer 
independent of p c , but obviously we recover the the New- 
tonian value in the Newtonian limit (p c — > 0, point B 
in Figures || and ||). Note that in general relativity all 
static n — 3 polytropes are unstable to radial perturba- 
tions (dM/dpc < 0). 

For each central density p c , there is a unique (maxi- 
mum) angular momentum, at which the equator reaches 
the Kepler frequency, and at which mass-shedding sets in. 
We mark this mass-shedding curve by the short dashed 
line in Figures || and ||. 

In order to find a transition from stable equilibrium 
to unstable equilibrium, we need to construct sequences 
of constant angular momentum and locate their turning 
points (Friedman, Ipser & Sorkin, 1988). We have marked 
sequences of constant angular momentum with thin, solid 
lines in Figures || and 0. They all connect the mass- 
shedding sequence with point B. The later can be under- 
stood by observing that, from equation d36j), T/j_W| scales 
with R^ 1 for constant j. As p c — > 0, the mass M remains 
finite, and hence R — > oo. We therefore have T/|W| — » 0, 
so that rotation plays an increasingly negligible role and 
we recover, for any j, the static limit. 

The sequences of constant angular momentum have a 
turning point only for a small range of angular momenta 
(18.25 < J < 20). Above these values, the sequences are 
monotonically increasing, while below they are monoton- 
ically decreasing. The thick solid line in Figures | and | 
connects the turning points, and therefore separates a re- 
gion of stable equilibrium (above the line), from a region 
of unstable equilibrium (below the line). 

As it cools, a supermassive star evolves toward higher 
density along the mass-shedding sequence and is therefore 
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Table 2 
Critical values at point A. 





Pc 


M 


Rp/M 


J/M 2 


T/\W\ 


Rp/ R e 


^/^spher 


1 T I T b 
1 j - 1 sprier 


Analytical Model 


6.6 x 1CT 9 


4.5551 


456 


0.8757 


7.44 x 10~ 3 


0.6667 


0.5443 


i 


Numerical Value 


7 x 1CT 9 


4.57 


427 


0.97 


8.99 x 10~ 3 


0.664 


0.5441 


1.15 



-n sphm . = (M/i?3)V2 

b / sphcr = 1 0.1130 R 2 p M 




Fig. 1. — Contours of constant mass energy density p in a plane containing the axis of rotation for a n = 3 polytrope at the critical point 
A (see Figs. |l] and ^). Each contour marks a change of the density by a factor of lO 1 / 3 . 



Mass -shedding Sequence 




Static Limit (TOV) 



Fig. 2. — Mass versus central density plot for relativistic, rotating n = 3 polytropes. The long dashed curve is the TOV solution for non- 
rotating, static configurations, and the short dashed curve marks the mass-shedding limit. The thin solid lines denote sequences of constant 
angular momentum, ranging from J = 15 (lowest curve) to J = 24 (highest curve) in increments o£ A J = 1. Turning points of these curves 
mark the onset of instability. The thick solid line connects these turning points (see also Figure H) and hence separates a region of stable 
configurations (above this line) from a region of unstable configurations (below this line). In particular, all nonrotating n = 3 polytropes are 
unstable to radial perturbations. A configuration evolving along the mass-shedding sequence with increasing central density becomes unstable 
at the critical point A. All sequences of constant angular momentum ^onnect the mass-shedding limit with point B for p c —f-Q (and hence 
_R — > oo). The mass of this configuration should agree with the mass M = 4.5525 of a Newtonian n = 3 polytrope (equation (|26|)), which we 
have marked by the open circle. The deviation of the solid point B from the analytical value is a measure of our numerical accuracy. 



4.57 



tS 4.56 - 



4 55 - 




Fig. 3.— Blowup of Figure ti. Here, the thin solid lines denote sequences of constant angular momentum raging from J = 17.5 (lowest 
curve) to J = 21 (highest curve] in increments of A J = 0.25. 
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stabilized by rotation until it reaches point A. At this 
point, the star becomes unstable to radial perturbations 
and will collapse. We summarize the numerical param- 
eters characterizing point A in Table 2, where they are 
compared with t he fi ndings of the analytical model calcu- 
lation of Section 3.2. Figure [j] shows the density profile 
at the critical point A in a plane containing the axis of 
rotation. 

Note that the numerical accuracy of an individual stellar 
model is fairly high. From the agreement of the limiting 
value B (solid dot in Figures || and ||) with the theoretical 
value (open circle), we estimate the numerical errors to be 
much less than 1 %. Identifying the critical point A, how- 
ever, requires locating the turning points along fairly flat 
curves, which therefore introduces a larger error. Never- 
theless, we expect that we can identify the characteristical 
parameters of point A to, typically, within less than 5%. 
Some of these parameters, however, change only very lit- 
tle along the mass-shedding sequence and can therefore by 
determined to much higher accuracy (for example T/|PF|). 
The number of digits in Table 2 reflects the uncertainty in 
each parameter. 

Having constructed self-consistent models of rotating 
n = 3 polytropes, we can now evaluate the quality o f th e 
predictions of the Roche approximation (see Section |3.l|) . 
We find that the predictions concerning the envelope alone, 
namely the ratio between equatorial and polar radius and 
the orbital frequency at mass shedding, are very accurate 
(Table 2). Figure [j] also shows that, at mass-shedding, only 
the envelope is distorted, while the core of the star remains 
nearly spherical. To predict T/\W\ at mass-shedding, 
the Roche approximation assumes the moment of inertia 
as well as the potential energy of the distorted star to 
be that of the corresponding spherical star, which yields 
T/|W| = 0.00744, independent of p c . In reality, the mo- 
ment of inertia will be larger for the rotating star (by about 
15% for configuration A, see the last column in Table 2), 
and the magnitude of the potential energy slightly smaller. 
Accordingly, we find that the Roche approximation under- 
estimates T/\W\ by about 20%. We also find that T/\W\ 
is not independent of p c (see Figure^). However, over our 
region of interest, T/|W| changes by much less than 1%, so 
that it can very well be approximated by 2 T/|W| = 0.009 
(cf. Table 1 in Hachisu 1986). 



3.4. Numerical Values in cgs units 

To obtain numerical values in cgs units from the dimcn- 
sionless "barred" quantities, we have to restore the poly- 
tropic constant K as well as c and G. Using eq. (||) for K, 
we find for the density p 



p c = K 3 p c 



= 9 x 10 10 gcm- 3 



(39) 



M 



7.0 x 10" 9 



where 7.0 x 10~ 9 is the approximate value of the central 
density of the critical configuration A. The temperature 
can now be found, for example, from eq. (17.2.9) in Shapiro 

2 Note that the relativistic definitions used to calculate T and W 
and (65) 



& Teukolsky (1983) 



T r = 9 x 10 iu K 



M 



-1/2 



Pc 



7.0 x 10" 9 



1/3 



(40) 



For SMSs with masses M > 10 6 A/©, the central temper- 
ature is always less than 6 x 10 7 A'. According to Fowler 
(1966), this is the minimum temperature required for gen- 
erating the SMS's luminosity via the CNO cycle. This 
justifies our assumption that nuclear reactions can be ne- 
glected. 

4. QUASISTATIONARY EVOLUTION TO THE ONSET OF 
INSTABILITY 

In this Section we discuss the evolution of a rotating 
SMS up to the onset of instability. In particular, we derive 
analytic expressions for the evolution of the mass, radius 
and angular momentum as a function of time. 

The evolution of the three quantities is not independent. 
Instead, they are coupled by two relations, namely the re- 
quirement that the star evolves along mass-shedding (so 
that T/|W| remains approximately constant), and that 
the angular momentum loss can be computed from the 
mass that leaves that star from the equator with the crit- 
ical angular velocity. This means that we can express, for 
instance, the angular momentum loss J and the change of 
radius R in terms of the mass loss M . The only relation 
that is yet to be determined fixes M and hence the overall 
timescale for the evolution of the three quantities. 

In a similar calculation, Bisnovatyi-Kogan, Zel'dovich & 
Novikov (1967) estimated M by constructing an approxi- 
mate stellar wind model, which they joined onto the outer 
envelope of the star. Their model depends on several un- 
known nondimensional parameters dealing with the wind 
solution. Here, we take a much more naive approach, de- 
termining the mass loss rate from the requirement that the 
star remain in equilibrium as it cools in a quasistationary 
manner. According to eq. (Q), the mass loss rate is related 
to the change of K (or the entropy) . and is thus governed 
by the star's luminosity (see eq. (E5J) below). 

As we have seen in Section 3.2, the post- Newtonian cor- 
rections and rotational contributions to the energy func- 
tional ( p0| ) are important for determining the stability of 
SMSs, but they have a very small effect on the equilibrium 
structure and can therefore be neglected for the purposes 
of this Section. Accordingly, the mass M of the star is well 
approximated by the Newtonian expression (Q) , which only 
depends on the polytropic constant K. The time deriva- 
tive of the mass is given by 



M _ 3 K 
M ~ 2 K' 



(41) 



The change of the total entropy S of the star is related the 
luminosity by 

L = -TS. (42) 

Using the first law of thermodynamics, the right hand side 
can be rewritten 

TS = J {Ts)dm = K J np 1/n dm = U—, (43) 

numerically are given in Cook, Shapiro & Teukolsky (1992), eq. (56) 
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Fig. 4. — T/|W| as a function of central density along the mass-shedding line. The dot marks the critical configuration A (see Figs. ^and^). 
The numerical calculation yields a value of T/[W] which is different from the value (T/|W|) shedd = 0.00744 that the Roche approximation 
predicts. Note, however, that T/|W| changes by far less than a percent over the region of interest, and can therefore well be approximated to 
be constant. 



where U is the internal energy (see Baumgarte & Shapiro 
1998, eq. [A38]). For n = 3 polytropes in equilibrium, the 
internal energy is equal to the negative of the gravitational 
potential energy 



2 R 



(44) 



(see eq. (^l"l), when post-Newtonian and rotational correc- 
tions are neglected). Collecting terms, we now find 



M __ L R 
M ~ ~M M' 

The luminosity L can be written as 

L = AivEdd, 



(45) 



(46) 



where the Eddington luminosity Z/Edd has been defined 
in As we have shown in Paper I, the luminosity of a 
SMS, rotating at the mass-shedding limit, is reduced by 
about 36% 

Ashcdd = 0.639. (47) 



Inserting (|4q) into (45), we now find 



M __ 4ttA R 
M ~ «T M' 



(48) 



As we have seen in Section [|, a SMS that evolves along 
the mass-shedding limit conserves T/IV4 7 ! to very good ap- 
proximation: 

T 

— ~ const. (49) 



Since T/|W| is proportional to J 2 M 3 R 1 , this implies 

(50) 



— = 2 3—. 

R J M 



Assuming that the mass-shedding occurs at the equator, 
we can write 



J = l csc M, 



(51) 



where ^ csc is the specific angular momentum of the escap- 
ing matter 



^csc ^shcdd-^ e . ^shedd-^p- 



(52) 



Since J can be written 



we find 



= \ 0.1130 R 2 p MQ she<id , 



J /3\ 3 1 M 



J (6 
1 ~ \2 



0.1130 M 



Inserting this into (KG) yields 



(53) 
(54) 

(55) 



Eqs. ([18|), ( p4D and ( p5| ) now completely determine the 
evolution of M, J and R. The equations can be simplified 
by defining the coefficients 



R M \\2J 0.1130 



k M = 47rA shcdd = 8.03 

/ 3\ 3 1 

kj=(-) = 29.86, 

J \2J 0.1130 

k R ee 2k, - 3 = 56.73, 



(56) 



in terms of which eqs. Jig), ( p4[ ) and ( p5| ) can then be 
rewritten 



(57) 
(58) 



M 


k M 


R 


M 


K 


M 


R 
R 


— ko — 

H M 




j 






1 


= kj M 





(59) 



The last two equations can be integrated immediately 



R 



Rc 



M 



J ( M 



Jcrit \ -^crit y 

Inserting ( |60| ) into ((5?]) then yields 

khlRcmt ( M 



M = -- 



K V ^crit 



M 

M crtt 



M" fc W = -^^ / ,11. 



(60) 
(61) 

(62) 
(63) 
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which can be integrated to 
M (. (k R - l)k M 



= 1 



l/(l-fcfl) 



(64) 



where — oo < t < and t = corresponds to arrival at the 
critical configuration. Defining 



tr 



(k R - 1) k M 
the last equation can be rewritten 



(!) 



M 



1 



t 



tr 



l/(l-fcfi) 



(65) 



(66) 



Note that i cr it is independent of the mass of the star, and 
determines the evolutionary timescale of SMSs. In cgs 
units, it takes the value 



(GM\ 



G (k R — 1) fcjVf ^ C 2 R ) crit 



8.8 x 10 n s 



(67) 



(compare Shapiro & Teukolsky 1983, where a similar 
timescale has been derived for SMSs that are stabilized 
by gas pressure rather than rotation). 

Similar expressions for R and J can be found by insert- 
ing (||) into @ and @: 



R 

Rent 

J 
J crit 



t 

^-crit 

t 



k R /(l-k R ) 



kj/(l-k R ) 



(68) 
(69) 



Using the identity k R = 2kj — 3, we find that the dimen- 
sionless quantities J/M 2 and R/M evolve according to 



J/M 2 



(J/M 2 



= 1 



Jcrit 



R/M 



= 1 



t 

^crit 

t 

tc.T\t, 



-1/2 



(70) 
(71) 



(R/M) cHt 

We plot M , J/M 2 and R/M as a function of time in Fig- 
ure H Note that J/M 2 and R/M evolve on the timescale 
icrit) which is why we adopt this value for our estimates 
in Section ^. The mass M and similarly K, however, ef- 
fectively evolve on a much longer timescale, which is a 
consequence of the small numerical value of the exponent 
1/(1 - k R ) = -0.0179. 

5. THE OUTCOME OF COLLAPSE: SPECULATIONS 



of the collapse of SMSs have been performed for nonrotat- 
ing configurations in spherical symmetry (see, for example, 
Shapiro & Teukolsky, 1979), for which black hole forma- 
tion is well established. These calculations, however, do 
not shed much light on the problem at hand, since rotation 
may play a crucial role in the dynamical evolution. The 
issue of the final fate of a collapsing, rapidly rotating SMS 
can therefore only be resolved by nonspherical, relativistic 
numerical simulations. 

In the meantime, however, we can attempt to assess 
crudely whether this collapse can actually lead to the for- 
mation of a supermassive black hole. To do so, we assume 
that each mass shell conserves its angular momentum dur- 
ing the collapse (cf. the discussion at the end of Section |J) , 
and consider two criteria. 

First, a particle can only be captured by a black hole 
if it is not be repelled by the angular momentum barrier. 
For simplicity, we take the newly formed black hole to be 
a Schwarzschild black hole, in which case a particle is in a 
capture orbit if 

I < 4m, (72) 

where / is the particle's specific angular momentum, and 
m the mass of the black hole. Here we ignore the effects of 
pressure, assuming that in relativistic collapse the matter 
approaches the nascent black hole supersonically by the 
time it enters the strong field domain. 

Second, for any portion of the star to form a Kerr black 
hole, the angular momentum J of that portion cannot be 
larger than the square of its mass m 



J/m 2 < 1. 



(73) 



We crudely evaluate these two criteria neglecting post- 
Newtonian corrections and nonspherical distortions due to 
rotation 3 . 

The specific angular momentum of a particle at radius 
r in the equatorial plane, at the onset of instability, is 

I = r 2 fUcdd, (74) 
where the critical orbital velocity Q cr n can be written 

°— =(a) bO =(3) (mL^- (75) 

Here, (i?/M) crit ~ 450 is the uniquely determined value 
at the onset of instability (see Table 2). Inserting the last 
two equations into (|72|), we find 



m(r) 1 /2\3/2/i?y/2M 
r 2 ~ I V3 ) VM Jcrit R 2 '' 



(76) 



where we have assumed that the mass of the black hole, 
rn, is the mass enclosed by a sphere of radius r. 

It is convenient to express the density and radius in 
terms of the dimensionless Lane-Emden functions 8 and £ 



Pc 



Once the supermassive star becomes unstable, it will 
start to collapse on a dynamical timescale. The outcome 
of this collapse depends on many factors, and can be de- 
termined only with a numerical, three-dimensional hydro- 
dynamics simulation in general relativity. Such a calcula- 
tion is only now underway (Baumgarte, Shapiro & Shibata 
1999). The only dynamical, fully relativistic simulations 

3 In reality, interior distortions and angular momentum exchange may play an important role late in the collapse, in particular if bars form, 
as we argue below. 



where a has units of length and is defined by 
( (n + Y)K P r n)/n \ 1/2 



(77) 



(78) 
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Fig. 5. — Evolutionary tracks of the mass M (solid line), angular momentum J/M 2 (short dashed line) and radius R/M (lone 
all normalized to their critical values at the onset of radial collapse (configuration A). The critical time t CI n is defined in eq. (f 
corresponds to the time at which the star reaches this onset of instability. 



dashed line), 
3), and t = 




~ Hiii i I i I i I i b 

0.2 0.4 0.6 0.6 1 



Fig. 6. — Plots of the derivative of the Lane-Emden function #(§), the enclosed mass rra(£) and the enclosed angular momentum J(£) divided 
by m(£) 2 as a function of £ for n = 3. The horizontal, dotted line in the top panel marks the threshold for capture by a central black hole, 
and the dotted line in the bottom panel marks the maximum angular momentum parameter for a Kerr black hole. 



(see, e.g., Shapiro & Teukolsky 1983). In terms of these 
quantities, the mass m(r) can be written 



m(r)=ATTa 3 p c e\0'(O\, 



(79) 



where the prime denotes a derivative with respect to £. 
For n — 3, the values of £ and 9' at the surface of the star 
are £i = 6.897 and |0'(£i)| = 0.0424. Using @ and ©, 
we find that (l76|) reduces to a criterion on #'(£) alone: 



no > v 2 f /2 r i? V /2 ^ 2 887 

-4V3J VMJcrit 2 - 887 ' 



(80) 



where we have adopted (i?/M) cr it ~ 450 for the critical 
configuration. 

We plot as a function of £ in the top panel of 

Figure 6, and mark the threshold value 2.887 by the hor- 
izontal line. A particle at a radius for which 6' (£) / '6' '(£1) 
is greater than this value can be captured, if the regions 
interior to this radius have collapsed to a black hole. For 
£ > £max ~ 0.57 £1, 0'(£)/0'(£i) is less than the capture 
threshold, and therefore the angular momentum barrier 
would prevent these regions from being caught by a newly 
formed, interior black hole. For the outermost regions, 
this is not surprising, since the configuration is critically 
rotating. Therefore, the outer region may remain in orbit, 
perhaps in a circumstellar disk, even if the rest of the star 
has collapsed. However, from the middle panel in Figure 6 
we find that this outer region only contains 5% of the 
mass, while about 95% of the mass could form a black 
hole. This, again, is a consequence of n = 3 polytropes 
being extremely centrally condensed. 

According to Figure 6, the angular momentum barrier 
would also prevent particles at very small radii, £ <J 0.05£i, 
from being captured by a black hole interior to that ra- 
dius. However, for this to be relevant, the initial black 
hole would have to be restricted to a very small fractional 



size, and as we will see below, the second criterion (73) 
does not allow such small black holes to form. 

We may reverse the above argument and view eq. ( |8C| ) as 
a condition on R/M. Since the left hand side, 0' (£) / 9' (£1) , 
has a maximum of about 6.6, R/M has to be smaller than 
about 2350 for black hole formation to be possible. It is 
interesting that R/M of the critical configuration A is only 
about five times smaller than this threshold compaction, 
just barely allowing the supermassive star to form a su- 
permassive black hole. This argument suggests that pri- 
mordial gas may have to pass through a phase as a SMS 
where it can lose angular momentum before it can possibly 
collapse to a SMBH. 

We can similarly evaluate the second criterion, eq. (f73"[), 
in terms of Lane-Emden functions. The angular momen- 
tum J(r) of the matter enclosed within radius r, rotating 
with the critical angular velocity, can be written 



J(r) = I(r) f2 s hedd 



(81) 



{ 2 s /2 R 1/2 M 8jL r 
v 3 ; y M' mt R 2 3 J F 



Dividing this by the square of the mass, eq. J7S|), then 
yields the condition 



1 > 



m(0 



^ M '^e\0'm 2 Jo 



In the bottom panel of Figure 6 we show the right hand 
side of this equation for a n — 3 polytrope and for 
(i?/M) cr it = 450. The dotted, horizontal line is the criti- 
cal value of unity. A region inside a radius £ can form a 
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Kerr black hole only if J(£)/m(£) 2 is less than unity. This 
condition is satisfied everywhere except for the innermost 
regions £ < £ m i n ~ 0.1 £1. At the surface, £ — £1, we 
recover the value (J/M) 2 rit = 0.876 found in eq. (|33|). 
Note that for any polytrope, m is proportional to r 3 for 



small radii r and J is proportional to m(r)Qr 



fir 5 . 



Therefore, J/m 2 scales like O/r close to the center, which 
prevents the formation of arbitrarily small black holes 
from rotating polytropes. In our case, the region inside 
£min ~ 0.1 £i, containing about 4% of the total mass, de- 
fines the minimum mass which can collapse to form an 
initial black hole. Since 9' (£) / '9' '(£i) exceeds the thresh- 
old for £ ;> £ m in, such a "minimal" seed black hole could 
accrete further material from the star. 

Finally, we note that a bar instability may form during 
the collapse of the star. The criterion for the formation of 
a bar on a dynamical timescale is 



T 

W\ 



£0.27 



(83) 



bar 



(see, e.g., Chandrasekhar 1969; Shapiro & Teukolsky 
1983). Since M and J are approximately conserved during 
the collapse, r/|W| scales with R~ l , so a bar will start to 
form when 



-Rbax = ( T \ (T_ 



0.009 



1 



Rent \ WV /.,, V |ir| / D.27 :!0 

The value of R/M at bar formation is then 

IV 
M 







fR\ 


^ bar 


30 1 





15. 



(84) 



(85) 



Eg. (g5|) suggests that the collapsing star may form a non- 
axisymmetric bar before it forms a black hole. This is 
an important result, since such a bar may result in a 
quasiperiodic emission of gravitational waves (cf. Smith, 
Houser & Centrella 1996). The frequency of these waves 
can be estimated from the expected bar rotation rate 
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For a SMS of 10 6 M Q this yields a frequency of 5 x 10~ 4 
Hz, which is in the range in which LISA is expected to 
be most sensitive (see, e.g., Thorne 1995). The frequency 
increases at later stages of the collapse, when M/R^ al be- 
comes larger. The strength of the gravitational wave signal 
can be crudely estimated to be 
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where Q is the star's quadrupole moment, and d is the dis- 
tance, which we scale to 1 Gpc (the Hubble distance is ~ 3 
Gpc). The signal strength increases with M/i?bar at late 
stages of the collapse. Apart from any bar, the implosion 
itself will be nonspherical, due to rotation, and will result 



in a gravitational wave burst. The collapse of a super- 
massive star may therefore be a very promising candidate 
for detection by space based gravitational wave detectors 
(LISA Pre-Phase A report 1995). 

Obviously, the arguments presented in this section are 
very crude and do not replace a fully self-consistent rela- 
tivistic hydrodynamical calculation. Our arguments nev- 
ertheless suggest that, upon reaching the onset of insta- 
bility, supermassive stars may form a supermassive black 
hole, containing a large part of the mass, and leaving only 
a few percent of the initial mass outside of the black hole, 
most likely in the form of a disk. We furthermore antici- 
pate that a bar may form during the collapse, which may 
lead to the emission of quasiperiodic gravitational waves. 

6. SUMMARY 

We have launched a fully relativistic study of the for- 
mation of SMBHs via the collapse of SMSs. In this paper 
we study the quasistationary, secular evolution of SMSs 
up to the critical configuration at which radial instabil- 
ity sets in and focus on the effects of rotation and general 
relativity. We identify the critical configuration and its 
characteristic parameters R/M, T/\W\ and J/M 2 . These 
parameters are independent of the mass of the SMS, and 
are therefore universal constants. The subsequent implo- 
sion, starting from this universal critical configuration, is 
therefore also uniquely determined and should produce a 
unique gravitational waveform. We compare results from 
an analytic, approximate treatment and a fully relativistic, 
numerical calculation, and find good agreement. We fur- 
thermore solve analytically for the time evolution of these 
parameters up to the onset of instability. 

Identifying the critical configuration at the onset of in- 
stability is interesting for its own sake. More importantly, 
however, this configuration will be adopted as initial data 
for future numerical simulations of SMS collapse (Baum- 
garte, Shapiro & Shibata 1999). In this paper we assemble 
qualitative arguments to anticipate the outcome of the col- 
lapse and find that the formation of a SMBH containing 
an appreciable fraction of the mass is not ruled out. Our 
arguments suggest that a transient phase as a SMS may 
be an efficient way for primordial gas to lose sufficient an- 
gular momentum in order to overcome the Kerr angular 
momentum barrier to forming a SMBH. 

Loeb & Rasio (1994) have emphasized that SMBHs ap- 
pear to have a minimum "observed" mass of ~ 10 6 M Q . 
Our crude arguments may help to explain this minimum 
in the context of the SMS formation scenario. SMSs with 
masses less than this value will not evolve in a quasista- 
tionary manner up to our critical configuration A, because 
some of our assumptions break down. For example, nu- 
clear burning may cause the star to explode (Appenzeller 
& Fricke 1972, Fuller, Woosley & Weaver, 1986). Electron- 
positron pairs may also be important for small masses, 
and may destabilize the star (see Zel'dovich & Novikov 
1971) before it reaches a small enough value of R/M to 
overcome the angular momentum barrier. On the other 
hand, all SMSs with M > 10 6 Af Q will evolve in a quasis- 
tationary manner until reaching the critical configuration 
A. The subsequent collapse will likely give rise to a SMBH 
of almost the same mass as the progenitor SMS. 

This paper was supported in part by NSF Grants AST 
96-18524 and PHY 99-02833 and NASA Grant NAG5-7152 
to the University of Illinois at Urbana-Champaign. 
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